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Abstract 

It is argued that colour octet contribution in exclusive charmonium decays is very impor- 
tant even though the same infrared divergence found in P-wave colour singlet inclusive 
decay is absent. This comes about because of the suppression at the level of the wave- 
function of P-wave quarkonia. Using the modified hard scattering approach of Brodsky, 
Lepage, Botts, Li and Sterman, we use charmonium decay into baryon-antibaryon pair as 
an example and show that the colour singlet contribution alone is clearly insufficient to 
explain the experimental decay widths. 

1 Introduction 

It is known that in inclusive decay of P-wave charmonium into light hadrons, there is 
an infrared (IR) divergence 1 in the reaction cc — ► gqq which is the leading inclusive 
contribution for xi but next-to- leading for xo an d Xi- This IR divergence reveals itself in 
the limit that the gluon becomes soft and the heavy cc are allowed to go on the mass-shell. 
In view of the divergence and the cc is a bounded system, one is obliged to keep them 
off-shell by an amount of the binding energy e. The resulting inclusive decay width has 
a logarithmic dependence on this binding energy r oc hie. However, it was shown 
in Ref. 2 that, in fact, this infrared divergence is there only due to the neglect of the 
next higher Fock state of the P-wave charmonium, the so-called colour octet component, 
which becomes degenerate with the valence cc singlet state in the already mentioned decay 
process in the soft gluon limit. 

In exclusive P-wave charmonium decay, this same IR divergence is not present because 
it only happens in a small region in momentum space and the decay products are again 
bound states. Therefore a priori, there is no good reason to introduce the colour octet. 
Indeed, the exclusive decay widths for xj have been calculated for a number of decay 
modes using the colour singlet valence Fock state alone and claims have been made that 
these calculations can account for the corresponding experimental partial widths. In 
the following, using the example of xj — * NN, we shall argue that, in fact, there is 
a very good reason to include the colour octet component in exclusive decay even in 
the absence of the IR divergence as found in the calculation of the inclusive width and 
point out a fallacy of some of the previous works. We show that within a self-consistent 
perturbative scheme, the modified hard scattering approach, colour singlet contribution 
alone is insufficient to account for the experimental data. Thus, based on the theoretical 
argument to be presented, colour octet is required. 
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2 Colour Singlet vs Colour Octet 



The decay of the cc system through annihilation is a short distance process with an 
annihilation size, L, of the order of the inverse charmonium mass L ~ 1/M XJ << 1. So 
for S-wave decay, one needs essentially ip s {L ~ 0), the wavefunction at the origin. For a 
P-wave system, the wavefunction at the origin is zero, so instead one uses Lift' (L ~ 0), 
the derivative of the wavefunction at the origin weighed by the annihilation size. In 
momentum space, the contributions from the wavefunctions alone will provide a P-wave 
to S-wave ratio of kip p (k) / M XJ ip s (k) so P-wave is suppressed in comparison to S-wave by 
1/M XJ . This suppression is very important, as we shall see presently, for whether the 
colour octet component should be included. 

In the hard scattering picture of Brodsky and Lepage, the decay probability amplitude is 
given by a convolution of the distribution amplitudes and the hard scattering perturbative 
part 

M Xj ^ N n ~ ./W07V ® fN<j)N ® fxAx.1 ® T H ■ (1) 

The only dimensional quantities in the amplitude are the decay constants and some power 
of the large scale of the process in question, namely the charmonium mass, hidden in the 
hard part T#. All these must together make up the right dimension of the probability 
amplitude which has one mass dimension because it is related to the decay width through 

r-ji-W. (2) 

The decay constants in Eq. (|I]) all have mass dimension two for both colour singlet and 
octet contribution to the width because both the wavefunctions of the baryon and the 
colour octet component are three-particle wavefunctions while the colour singlet compo- 
nent is a P-wave two-particle wavefunction. Using the right power of the charmonium 
mass from the hard part T# to make up for the right mass dimension of the amplitude 
gives 

As one can see, both the octet and singlet contribution depend on the same power of 
M XJ . The octet contribution is not negligible as far as the large scale M XJ is concerned. 
This is the case only because, as we saw above, P-wave charmonium is weighed down by 
the annihilation size L ~ 1/M XJ in the amplitude. Had we been dealing with the S-wave 
charmonium J/ip instead, the singlet amplitude would go as and the colour octet 

can be neglected in comparison with the colour singlet valence component. Returning 
to the xj decay, possible suppression from the a s dependence is absence. The singlet 
and octet contribution are only differed by ^foT s for typical graphs which is not small. 
Therefore there is no good reason to neglect the colour octet from exclusive decay. 
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3 Colour Singlet Contribution within the Modified 
Hard Scattering Approach 

In the modified hard scattering picture, the decay probability amplitude is given by 

M xj^nn ~ ipB ® ifis ® ip X j ® t h® exp(-S) . (5) 

a convolution of the wavef unctions, hard part and Sudakov factor. In the valence colour 
singlet contribution, the part of the convolution that can potentially render the result 
ambiguous is the nucleon wavefunction. The reason being that there are many nucleon 
model wavefunctions available that one can use, for example see Ref. 3, for this calcu- 
lation. However, none of them can explain data on nucleon form factor at low Q 2 < 50 
GeV 2 , which is the relevant energy range of interest in our present problem. As shown 
in Ref. 4, these models yield nucleon form factor results that are well below the experi- 
mental data. As a result, one cannot possibly hope to use these models to calculate the 
Xj partial decay width. In view of the smallness of the perturbative contribution, one is 
forced to leave the perturbative region and find non-perturbative contribution which can 
account for the data. The full expression for the nucleon form factor can be written as 

oo „ n 

Fi(Q 2 ) = E J H[dV] E e i ^ n \x, k±) *4) (6) 

n=3 je{q, q} 

where k' ± = k± + (1 — Xi)q for the struck quark and k' ± = k± — Xjq for the spectator 
quarks, the second sum is over the parton index of the quark and antiquark constituents 
and we have only put a representative x and kj_ to stand for all the momentum fractions 
and internal momenta of the constituents in the argument of the n-constituent nucleon 
wavefunction ip^ n \ When hit by a highly virtual photon, one knows that the perturbative 
contribution comes from the region where the internal momenta of either one of the 
wavefunction in Eq. (^) are hard while those in the other remain soft. Since this is not 
the dominant contribution, one looks at the region where the internal momenta in both 
wavefunctions are soft. The reason being that bounded hadrons tend to have small internal 
momentum so if the incoming photon momentum is not so hard, as in the low Q 2 but 
Q 2 » MeV 2 region, the overlap of both wavefunctions with small internal momentum 
should dominate. This can happen only if the struck quark takes up almost all the 
momentum of the nucleon so that (1 — xi)q and Xjq all become soft. This is most unlikely 
if many constituent partons are there to share the nucleon momentum. So one can take 
only the valence Fock state as the dominant contribution 




This simplication enables one to fit this to the F\ data, GRV valence quark distribution 
and other data like r(J/ip — ► NN) to obtain a model wavefunction for the nucleon 5 
valid in the low Q 2 range of interest. The other remaining colour singlet xj wavefunction 
is more or less standard (see for example Ref. 6). The hard part T# in the MHSP comes 
from 4 graphs up to permutation of the light quark-antiquark lines 7 . Here we do not give 
the full expression but just to point out that in the MHSP, the probability amplitude has 
the following dependence on a s 

~ a s (ti) a s {t 2 ) a s {t 3 ) (8) 
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J 


rW(Xj-flP) [eV] 


PDG 8 [eV] 


1 


7.49 


75.68 


2 


46.97 


200.00 



Table 1: Colour singlet contributions to xj decay into pp. 

where the renormalization scales ij, % = 1,2,3 are determined by the largest scale in 
the neighbourhood of the vertices that carry the particular a s . The scales could be the 
virtualities of the neighbouring propagators or the inverse of the smallest transverse size 
squared between the quarks (antiquarks) in the baryon (antibaryon). So the renormaliza- 
tion scales are determined dynamically by the process which is a typical hallmark of the 
MHSP. The product of a s 's are therefore part of the integrand of the convolution. One 
cannot, as done in some earlier work, take this product outside the integrand and assign 
some constant effective value for the coupling. In our opinion, this procedure is rather 
unconvincing because A4^ ~ and therefore ~ oP s . Changing a s from a value of 0.3 
to 0.5 will make a factor of 20 difference in the width so provided the singlet contribution 
is of reasonable size, one can always choose a value such that the colour singlet alone can 
account for the experimental width. 

As we argued in the previous section, colour octet is not suppressed for P-wave decay 
with respect to the singlet contribution therefore both should be needed to account for 
the experimental data. In table p], we listed the singlet contributions to xj — ¥ PP partial 
width. Xo is forbidden in this decay mode by angular momentum conservation. One can 
see that the singlet contribution can only account for 10% of xi an d 24% of xi °f t ne 
experimental partial width. So in agreement with our theoretical argument presented in 
the previous section, colour octet is needed in exclusive P-wave charmonium decay even 
in the absence of the same infrared divergence as found in the colour singlet inclusive 
decay width, to cancel which the colour octet was introduced. 

Our theoretical argument and conclusion are however much more general than just for 
charmonium decay alone. One could apply these equally to bottomium as well as to P- 
wave quarkonium production. In these processes, colour octet should be required. One 
can however venture even further by making the suggestion that whenever a P-wave heavy 
quark- ant iquark pair, which could be a part of the constituents of some component of a 
hadronic wavefunction, is needed, one will also have the need of the component of the 
corresponding colour octet accompanying the same set of other constituents as well. 
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